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We show that the Calogero-Sutherland (C-S) gas, a famous exact soluble one-dimensional system 
with an inverse square long range interaction, can be realized by dimension reduction in a cold 
Bose atom system with a dipole-dipole interaction. Depending on the orientation of the dipoles, 
the effective interaction is either attractive or repulsive. The low-lying effective theory may be a 
Luttinger liquid when the exclusion statistics parameter A may be well-defined. We hope that the 
C-S gas can be realized experimentally and the Luttinger liquid character can be observed. 
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Introduction One-dimensional ultracold atom systems 
have shown amazing physical behaviors in a bunch of 
experiments in recent years. A phase transition from 
a one-dimensional Bose Einstein condensate(BEC) to a 
Mott insulator has been observed in a one-dimensional 
optical lattice Q . The strong coupling limit of the one- 
dimensional gas with a contact interaction, known as the 
Tonks-Girardeau gas , has also reached in the ultracold 
Bose atoms 0- 

Another kind of important one-dimensional systems 
with inverse square interaction was known as the 
Calogero- Sutherland (C-S) models HH]. Nowadays the 
C-S models play a fundamental role in the contemporary 
theoretical physics. The eigen wave functions of these ex- 
actly soluble models may be explicitly expressed by sym- 
metric polynomials [6|, |3 ■ The quasi-particle excitations 
in these models obey exclusion statistics 0, 0, liol. The 
low-lying effective theory is a Luttinger hauid [ill 
characterized by the exclusion statistics parameter A 
0, 0| . It was shown 0, Q| that this Luttinger liquid 
is equivalent to a c = 1 conformal field theory on a circle 
with a radius R = •\/2/A 0|. The C-S model may be a 
universal Hamiltonian of some weakly disordered metals 
[l^. The models may also be used to describe the edge 
excitations in the fractional quantum Hall effect [l^ . 

Besides mentioned above, there were already numer- 
ous works in studying the C-S type models [Til' , but un- 
til now a one-dimensional physical system with a real 
inverse square interaction still lacks experimental real- 
ization. Recently, a major progress in the cold atoms 
was that chromium atoms ^^Cr with a dipole-dipole in- 
teraction were turned into a BEC j3| ■ This provides an 
opportunity to induce an inverse square interaction by 
the dimension reduction. There already are many theo- 
retical considerations for the BEC with dipole-dipole in- 
teractions in ultracold gases [23,0|. To our knowledge, 
such a dimension reduction to the dipole interaction was 
not studied, which will be a main topic of this Letter. 
Reducing to a single-hand model With a very strong z- 
direction confinement, a cold atom BEC cloud has a 
pancake-like shape trapped in the x-y plane. Using an 
periodic optical potential along the x-direction, the pan- 



cake is incised into an array of cigar-like BECs. Apply 
an external field in (or perpendicular to) this x-y plane 
so that the dipoles in (perpendicular to) this plane are 
parallel. (See Fig. 1(b) or (c). We first neglect the slow 
varying harmonic trap in the a;-direction.) This is a prob- 
lem of the Bose atoms loaded into an optical lattice with 
high filling factor. For the alkali metal atoms, the sys- 
tem is described by multi-band Bose Hubbard model but 
can be mapped into an effective single-band model [j^ . 
If the product of the average number A^o in a BEC and 
the on-site repulsion U is much larger than the trapping 
potential along the y-direction while it is much smaller 
than the strength of the optical lattice, the wave func- 
tion may be approximated by the product of the single 
atom ground state wave function in the x-direction and 
the BEC's in the y direction. In this approximation, the 
atomic field operator may be written as 

ipir) Ki^aiw(x - Xi)'^TF{y) (1) 

i 

where w{x) is the Wannicr function along the x-direction 
and ^TF{y) is the Thomas- Fermi wave function along 
the 2/-direction, which describes the ground state of 
atoms in each single site. a\ (a^) is the creation (an- 
nihilation) operator for an atom located in the ground 
state at each site. The rcnormalized hopping tn and 
interaction Uji can be estimated in the Thomas- Fermi 
approximation 22] . The hopping t may not be rcnormal- 
ized while Ur is reduced from its bare value con- 
siderably due to the repulsive on-site interactions that 
spread out the condensate wave function. The consis- 
tent condition for this reduction of Ur is Ip/a <^ Nq <C 
{hujy/'hujp)'^\/2TTLy/a, where a is the s-wave scattering 
length, hojy is the trapping energy in the j/-direction, 
Ly — ^JKJmWy is the effective length in the y-direction. 
Ip is the lattice spacing and hio^ = Ti' jml'^. These condi- 
tions may easily reach in experiment. 

For the atom cloud with dipole-dipole interaction, be- 
sides the on-site interaction, there is another interaction 
between two parallel dipoles located in r and r' , which 
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reads 



(2) 



where 9 is the angle between R = r — r' and d is the 
magnitude of the dipole moment d . For dipoles laying 
in the plane, Q is shown in Fig. 1(a) while Q — Tr/2 
for dipoles perpendicular to the plane. The Thomas- 
Fermi density of dipolar BEC has the same parabolic 
form as that in the short range interaction BEC. The 
dipole-dipole interaction only affects the central den- 
sity and the Thomas- Fermi radius Therefore, the 
wave function factorization approximation is also valid if 
hujy < NoU[l + 0{edd)] < hwp where U[l + 0{edd)] is the 
interacting energy scale li^. {edd < 1 for the practical 
system will be defined later.) 

Consider a system with atom number LNq+N{N < L) 
where L the x-direction lattice size. The effective single- 
band lattice model can be obtained by keeping only the 
nearest neighbor hopping and expanding on-site energy 
to the second order near mean occupation Nq jl^, HJ, 
i.e., 



H/N^ = - ^ ta\aj + ^ ^ n,{n, - 1) 

(ij) » 

+ ^Ud,tjninj, (3) 



where Ui - 
Vp{x))w{x 
bor sites, 
and Ud 



No; Jdx nj*{x^x,){-^±,+ 
Xj) for a pair of nearest neigh- 

Ur = ^Jdx\wixrjdy\^TF{yr 

X^)\'^\w{x' - 



mNo 

:,y - W7jdrdr'\w{x 

x,)\^Vd{ry)\^TF{y)\'\^TF{yT. 

This is a homogeneous model by neglecting the trap- 
ping in the a;-direction and describes the dynamics of 
the fluctuation around the mean occupation A^o- Adding 
back the harmonic trap with the frequency ujx, the sys- 
tem has an effective size L — ^Jfij muj^ in the x-direction 
. If the width of a single BEC is wider, the number A^o 
of atom per site may vary site by site and Nq in the cen- 
ter of trapping may be several times of Nq in far away 
from the center. Within a cigar-like single dipolar BEC, 
the dipole-dipole repulsion is strong enough against ras- 
ing the width of the single BEC at the trapping center. 
Therefore, the factorization of the ground state wave is 
still valid and A'o keeps an constant. The resulting inho- 
mogeneous model is iJ in JSJ added by a term Vh^iJii 
where 



—mujz 
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dx\w(x — Xi)\ 
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X 



(4) 



is the harmonic trapping potential along the x-direction. 

We now assume the trap in the y-direction is very weak 
such that |\E'tf(2/)P can be approximated by the average 
density poi which is consistent with another approxima- 
tion we will use, i.e., assuming the condensate along the 



y-direction is infinitely long. Under these approxima- 
tions, it is easy to integrate over y and y' and arrives 
at 



Ud, 



Gxy{z) / dxdx'\w{x — Xi)\ \w{x' — Xj)Y 



G 



[x — .t')2 ' [xi — XjY 



(5) 



where Gxy — —2d^/5ocos^9o if the dipoles orient in the 
plane and Gz = Sd^po if dipoles is perpendicular to the 
plane. Here, Oq is the angle of the dipole with respect 
to the X axis ( FigQ^a)). If the on-site interaction is 
switched off, to prevent the BEC collapse due to the at- 
tractive dipole-diploe interaction in a single BEC, one 
takes cos^ > 5 such that the dipoles in the y-direction 
are repulsive. Thus, we have a one-dimensional lattice 
model with the on-site and inverse square interactions 
between the single modes from the ground state. If the 
dipoles lay in the x-y plane (Fig. 1(b)), the effective in- 
teraction between those atoms in different sites is attrac- 
tive while it is repulsive if the dipoles are perpendicular 
to the plane (Fig. 1(c)). 

The stability condition of a single BEC at a given site 
is given by edd = md^ / {2,Tt? a) < 1 For ^"^Ci, the 

experimental measured e^ti = 0.159±0.034 [23|. However, 
the renormahzed on-site interaction is reduced to l23| 



Airh a Ip 



(6) 



m Rtf 

where Rtf is the Thomas- Fermi radius of the single 
BEC. On the other hand, d is not renormalizable as 
shown in Thus, although edd < 1, the renormahzed 
on-site interaction may be much weaker than the dipole- 
diploe interaction since eddRTF /Ip 3> 1 as Ip/Rtf ^ 1 
for the cigar-like BEC we are considered. Therefore, 
in the effective one-dimensional single-band model, the 
dipole-dipole interaction dominates and the on-site in- 
teraction may be turned off. 

Dilute gas limit and CS model In the dilute gas limit with 
N/L <C 1, the A^-particle system may be described by a 
continuous model because the dispersion —t cos klp/Ti ~ 
k^/{2mji) with m^j ~ fi^/itlp). Neglecting the on-site 
interaction as argued above, the effective Hamiltonian 
reads 



CS 



N 

E 
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mRLj^xl 



,(7) 



where the renormahzed trapping potential is defined by 
TOW^ = mpuj'^. The model described by the Hamiltonian 
(fTI) is the famous C-S model 0, ■ A 'particle' at Xi 
in this continuous model corresponds to one more atom 
than iVo at the lattice site i in the lattice model. The 
ground state wave function is given by 

^oA^i,-,^n)= n (•^,-:^;fc)^e-*'™^^%(8) 

l<j<k<N 
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(a) 




(b) 





FIG. 1: (Color on-line) The trapped cold atom cloud in an 
optical lattice, (a) The parameters used in deriving eq. Q 
for the case (b). (b) The dipoles are in the x-y plane, (c) The 
dipoles are perpendicular to the x-y plane. The angle 9 is the 
same as shown in (a). 

where A is the solution of the equation A(A— 1) = —gxy = 
(mii/h^)Gxy or = gz — {mi^/h^)Gz, respectively. The 
ground state energy Eg = ^NfiuJii{X{N — 1) + 1). 

For A = A:rjy = (1 - y/l - 45^y)/2, if g^y < 1/4 the 
ground state © is well-defined because Xxy > is real . 
Although the interaction is attractive, the particles in one 
dimension mutually exclude because the Jastraw factor 
in the ground state wave function forbids particles occupy 
the same position. The ground state energy Eg will be 
imaginary ii gxy > j- This means that the ground state 
is not stable when gxy > j ■ The collapse of the ground 
state reflects the factorization of the wave function is false 
if the attraction between atoms is too strong. The system 
is back to a two-dimensional BEC if edd < 1 because the 
on-site repulsive is not renormalized. 

For the repulsive coupling constant gz, X± = 
i±Vi+4g. ^ Pqj. _ 1 < a_ < 0, i.e, gz <j, the wave func- 
tion is still square integrable. This implies if the repul- 
sive between particles is not strong enough, the particles 
would like to gather together in the real space. Hence, 
the long wave length approximation used here is not valid 
and one has to be studied in the lattice model. We do 
not go this complicated case here. For a strong repul- 
sion with 5z > f , ® is no longer square integrable for 
A = A_ < — |. Instead, one should take A = A+ > 1.5 in 
This is the standard solution of the C-S model. 

From now on we will focus on the C-S gas, i.e., gxy < \ 
and gz > J- For the C-S gas, all the excitation states may 
be constructed and may be expressed by so called hidden- 
Jack polynomials '2^. We can, in principle, calculate any 
physical observable. However, to directly measure these 
observable in a trapped gas is difficult. In experiments, 
a well-established technique is switching off all trapping 
potentials and letting the atom cloud freely expand. Af- 



ter a long enough time, record the image of the cloud. 
In this time of flight, the image of free atom cloud re- 
flects the momentum distribution of the trapped atoms. 
Due to the harmonic trap, the momentum is no longer a 
good quantum number. Thus, the calculation from these 
exact eigen states in the thermodynamic limit becomes 
unwieldy. A frequently used approximation is the local 
density approximation. The starting point of this ap- 
proximation is solving the homogeneous system in the 
thermodynamic limit, i.e., lor = and N and L 00 
while p = N/L is fixed. Solving the model in this limit 
yields to solve the thermodynamic limit of the Sutherland 
model with a period potential 



V ^ 



■E 



7r2A(A - 1) 



2mR ^ L^sin [TT{xi - Xj)/L] 



(9) 



The quantum states of the Sutherland model are la- 
belled by the pseudo-momenta which are determined 
by the Bethe ansatz equations Lki = 27r/i -|- 7r(l — 
•^) X]j<i sgn(A;j — ki). Due to the exclusion between the 
particles, the ground state is a pseudo- Fermi sea which is 
given by {fcO} = {-^(iV - l)-fiN - 3), fiN - 
3), ■^{N — 1)} . The pseudo-Fermi momentum is kp — 

Notice that the situation in A < 1 is very different 
from A > 1. For A > 1, in a physical momentum in- 
terval there is only either one or zero one since the 
pseudo- momentum interval is 2jA > But for A < 1, 
as we have seen, the number of ones is larger than one in 
an interval Such a cluster behavior of the particles 
stems from the attraction between the particles if A < 1. 
To see this matter clearly, we calculate the single parti- 
cle correlation function. In fact, all correlation functions 
and their asymptotic forms were calculated. (See, e.g., 
0,^3-) the thermodynamic limit, we need only know 
the asymptotic ones. For the particle ^"^(2;) with an ex- 
clusion statistics parameter A, the single-particle correla- 
tion function in the a; — > cx) limit is given by G'(x, 0; A) = 

For A = A+ > 1, the mo- 



(VI.;,(X,0)VP1(0,0)) 

mentum distribution of the single-particle is well-defined 
and continuous near the Fermi momentum pp = hkp/X, 
G\+{p) = Gx+{pf) - const. • sgn(p - pp)\p - 
where p is the physical momentum with interval 
Thus, it is a standard Luttinger liquid [lllll^. 

However, for A = Xxy < 1, the momentum distribu- 
tion of the single-particle is divergent at the Fermi point 
pp, i.e., G\^^{p) sgn{p - pp)\p - ppl^'^y-^. This di- 
vergence implies that there is no single particle prop- 
agation. The instability of this single particle Fermi 
surface comes from the attraction between the parti- 
cles. The competition between this attraction and the 
zero-point motion may lead to a cluster propagating be- 
havior as we have analyzed in the pseudo-Fermi sea. 
In the coordinate space, this cluster behaviors can be 
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seen as follows. Defining a composite particle opera- 
tor <^nX.M = lim,^o[a-*"^"-^)*A.„(^)*A.„(x + 
a)...'^ \^y{x + {n — l)a)], for an integer n, the momentum 
distribution of this composite particle near is given by 
Gn\^y{p) = G„A,JpF)-const.-sgn(p-pi.)|p-pFr^"""^ 
The minimal n with nX^y — 1 > is n = [l/Aa;^] + 1. 
However, the composite particle with n' = [l/X^y] — 1 
will accompany with this n = [1/ Xxy] + 1 particle to in- 
sure the pseudo-Fermi sea is filled up but the former has 
a divergent distribution at the Fermi surface. A consis- 
tent choice IS n = I for Xxy = s/l and thus the exponent 
is IXxy — 1 = s — 1, e.g., Xxy = 2/5, s — 1 = 1. This 
reflects the exclusion statistics of the quasi-particles. 
Time of flight and local density approximation The time 
of flight experiment directly measure the momentum dis- 
tribution in the trapped gas. In the y-direction along the 
cigar-like BEC, the dynamics is controlled by Bogoliubov 
quasiparticles excitation while in the cc-direction, the dy- 
namics is dominated by the Luttinger liquid described 
above. We hope this Luttinger liquid behavior may be 
observed by measuring the slope of the density profile at 
the Fermi surface in the x-direction. The density-density 
correlation function {p{x, 0)p(0, 0)) ~ x~^ cor resp onds 
to the shot noise in the time of flight image |2y|. To 
experimentally observe this Luttinger liquid requires a 
very high resolving rate because there is background LNq 
atoms' image. 

For ujn 7^ 0, one may use the local density ap- 
proximation. Defining the local Fermi momentum 
Pf{x) — \/p>p — m?j^u)'j^x'^ / X"^ , the density distribution in 
the Thomas-Fcrmi approximation, which is valid in the 
thermodynamic limit, is defined by 27] 



I 



PF 



dP- 



X Gx{p,P)P 



A2 



PF 



Pi 



Gx{p,Pf) (10) 



where pf{L) = \/p% — mp fiujji/ X and is in between 
Pp and pp{L). Unlike the short range interaction Bose 
gas, the long range interaction gas in one-dimension is 
more like the degenerate Fermi gas: Despite the spa- 
tial anisotropy of the trap, the momentum distribution 
is uniform in the fca;-direction |27j . 

Conclusions We have shown the possibility to realize the 
C-S gas in a cold atom system with dipole-dipole inter- 
action. The ground state behaviors of the system in var- 
ious parameters were discussed. The experimental impli- 
cations to explore the properties of these ground states 
and low-lying excitations are suggested. 
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